Single-particle and Interaction Effects on the Cohesion and Transport and Magnetic 
Properties of Metal Nanowires at Finite Voltages 
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The single-particle and interaction effects on the cohesion, electronic transport, and some magnetic 
properties of metallic nanocylinders have been studied at finite voltages by using a generalized mean- 
field electron model. The electron-electron interactions are treated in the self-consistent Hartree 
approximation. Our results show the single-particle effect is dominant in the cohesive force, while 
the nonzero magnetoconductance and magnetotension coefficients are attributed to the interaction 
effect. Both single-particle and interaction effects are important to the differential conductance and 
magnetic susceptibility. 

PACS numbers: 61.46.+W, 68.65. La, 73.23.-b,85.70.-w 



I. INTRODUCTION 

Metal nanowires have been the subject of many exper- 
imental and theoretical studies^. An important feature is 
the quantization of motion of electrons because of spatial 
confinements. In the linear regime, the Coulomb inter- 
action among the electrons are not important, and the 
transport properties can be well described by Landauer 
formula— in the framework of the free electron model 
where the transmission probability can be calculated at 
equilibrium. Under a finite bias, however, the scattering 
states of right- and left-moving electrons in a nanowire 
are populated differently, even if there is no inelastic 
scattering within the nanowire. An adequate treatment 
of the electron-electron interactions is therefore crucial 
to correctly describe this non-equilibrium electron distri- 
bution in the contact. Some studies of cohesion^ and 
transpor1jii£i£iL& in metal nanowires at finite voltages us- 
ing continuum models did not include electron-electron 
interactions so that the calculated transport and ener- 
getics depended separately on both the left and right 
chemical potentials // + and thus violating the "gauge 
invariance" condition: The calculated physical quanti- 
ties should depend only on the voltage eV = n + — /Lt_, 
and should be invariant under a global shift of the bias 
since the total charge is conserved^ Other calculations 
have been made for non-equilibrium metallic contacts 
including the electron-electron interactions within the 
local-density approximationiiSiiiiiS However, these cal- 
culations utilized the canonical ensemble, which is not 
appropriate for an open mesoscopic system. Finally, a 
self-consistent formulation of transport and cohesion at 
finite bias has been developed based on ab initio and 
tight-binding calculationS ) 13 i 14 ' 15 i 16 ' 17 but ab initio cal- 
culations can thus far only simulate small-size systems. 

In this paper, we use the extended nanoscale mean- 
field electron model developed in Ref. to investigate 
the single-particle and Coulomb interaction effects on 
the cohesion, transport and magnetic properties of metal 
nanowires at finite voltages. Since those quantities are 



response functions, they can be used to characterize how 
the single-particle motion of electrons and their Coulomb 
interactions respond to the corresponding external forces. 
Those properties have been analyzed by Zagoskin^ and 
Bogachek et. am using a free-electron model which can 
only take into account the single-particle effect. Here 
we emphasize the necessity of an adequate treatment 
of electron-electron interactions not only to satisfy the 
gauge invariance, but also to show that the response of 
the Coulomb interactions to the applied voltage and mag- 
netic field can give rise important effects on the cohesion, 
transport and magnetic properties at finite voltages. 

This paper is organized as follows: In Sec. [HI we 
briefly introduce the extended nanoscale mean-field elec- 
tron model. The details of this model are described in 
Ref. ^3 The the cohesion, differential conductance, and 
magnetic properties are studied in sections IIIII and IIVI 
Sec. Ivl presents some discussions and conclusion. 



II. MODEL 

Here we briefly introduce the extended nanoscale 
mean-field electron model. Details can be found in Ref. 
ITsl In this paper, we only consider the temperature 
T = 0. We consider a cylindrical metallic mesoscopic 
conductor connected to two reservoirs with respective 
chemical potentials £t+(-) — fJ- + eU + (_), where [i is the 
electron chemical potential in the reservoirs at equilib- 
rium, and V+(_) is the voltage at the left (right) reservoir. 
While there is no general prescription for constructing a 
free energy for such a system out of equilibrium, it is 
possible to do so based on scattering theory if inelas- 
tic scattering can be neglected, i.e. if the length L of 
the wire satisfies L <C Lj n . In that case the scattering 
states within the wire populated by the left (right) reser- 
voir form a subsystem in equilibrium with that reservoir 
and the dissipation only takes place within the reservoirs 
for the outgoing electrons. Using the hard-wall bound- 
ary and treating the electron-electron interactions in the 
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Hartree approximation, one can define a non-equilibrium 
grand-canonical potential f2 of the system at zero tem- 
perature a o 1 ^ 1 ? 



Q [Ro,V,U] 



a=± 



3Ai 



3/2 



7V+C7, (1) 



where e F and k F are the Fermi energy and Fermi wave 
vector, respectively, e a = ji a — U, and e v = SF^lt/kpR^ 
is the v-th. transverse eigenvalue of an electron in a cylin- 
drical wire with radius Rq, with j u the zeros of Bcsscl 
functions, U is the Hartree potential energy, which must 
be self-consistently determined, and N+ is the total num- 
ber positive background charges and is taken to be 



N + = 



kpRgL 



kpRoL 



k F L 



(2) 
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The second term on the r.h.s of Eq. J2J) corresponds to the 
well know surface correction in the free-electron model^S 
which is essentially equivalent to placing the hard-wall 
boundary at a distance d — 2>i:/8k F outside the surface 
of the metal wireiSi The last term represent an integrated 
curvature contribution. 

Based on the assumptions of our model, the Hartree 
energy U is constant along the wire and can be self- 
consistently determined by the following charge neutral- 
ity condition at a given voltage V, 

Q = ^e[N-(fx+ -U) + JV_0i_ - CO] - eN+ = 0, (3) 



where 



2L 



{H± -U-e v ) 



1/2 



(4) 



is the number of right (left)-moving electrons in the cylin- 
drical wire up to energy £*+(-) — U. The summation is 
over all states such that e v < — U. The unifor- 

mity of the Hartree potential U can be destroyed either 
by the realistic atomic structure of wire (including im- 
purities in the wire), which can cause both elastic and 
inelastic scatterings, or the nonidcal couplings between 
the wire and the reservoirs, which induce backscattering 
(see Ref . Il8l for detailed discussion). For our purposes in 
this paper, these effects are not important. Equation © 
gives a relation^ 



U = U s + + - e F , 



(5) 



where U s is calculated with a symmetric voltage drop 
V + = —V- = \V between the two ends of the wire. 
Equation JSJ) will guarantee that all physical quantities 
calculated in the following are just a function of voltage 
V, and not of [i- and /i+ separately. 

The current in the cylindrical wire, according to our 
assumptions, is given as 



where h is the Planck constant. 
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FIG. 1: (Color line)The differential conductance and cohe- 
sive force of a metal cylindrical nanowire versus voltage for 
a nanowire with radius k F Ro = 10.64. Here G a and G u are 
the contributions from single-particle motion and interaction 
effect, respectively, and G = Gt is the total differential con- 
ductance given by Eq. 
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FIG. 2: (Color line) The derivative of the Hartree energy U 
with respect to the voltage. 



III. COHESION AND DIFFERENTIAL 
CONDUCTANCE 



The cohesive force F 



an[Ro,v,u] 

dL 



where the 



N 4 



derivative is taken at constant background positive 
charge N +1 is given by 



F(R Q ,V,U) = 



4£f \ - / £g - By 

3A F M ef 

a— ± x 
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£a - £ » + 3k F L- ll dRo 



e F 



k F R% dL 



(7) 



The first term comes from the single-particle levels and 
the second term is due to the finite size effect. The con- 
tribution from the derivative of Hartree potential with re- 
spect to L in the electronic free energy is completely can- 
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celed by that from the positive charge background, and 
the interactions affect the cohesive force only by shifting 
the chemical potentials, or equivalently by shifting the 
single-particle levels by an amount of the Hartree po- 
tential U . Therefore, the cohesion of metal nanowires is 
determined by the single-particle motion of the electrons, 
even at finite voltages. 

The differential conductance G 



|f is given by 



G 



G(R ,V,U) = ^J2 d ^- £ ^ 



a = ± 



dU 



d(eV) 



(8) 



where Go = 2e 2 /h is the unit quantum conductance. We 
see that both single-particle motion and interactions con- 
tribute to the differential conductance. 

The differential conductance G and cohesive force as a 
function of voltage V are shown in Fig. ^ In the figure, 
we have split the total differential conductance G into the 
single-particle contribution part G s and the interaction 
contribution part G u . At small voltages, G — G s and is 
equal to that obtained from free electron model, which 
is consistent with the result of linear transport. At large 
voltages, there are peak structures at the conductance 
jumps, which come from G u . This is because at the sub- 
band thresholds, a small voltage change can induce large 
fluctuation of the Hartree potential, as can be seen in Fig. 
121 From Fig. 2] one can also see the correlation between 
the variations of the conductance jumps and the force 
oscillation as a function of the voltage, which should be 
observable experimentally. 

We should mention that an adequate treatment the 
electron-electron interactions (Hartree approximation in 
this work), is essential to satisfying the gauge invariance. 
The calculated cohesive force and differential conduc- 
tance are just functions of the voltage between the two 
ends of the wire. This is in contrast to the results in Refs. 
and0, which depends on the partition (3 of the voltage 
drops between the two contacts by using a charge neu- 
trality: N-(£ F +l3eV) + N_(e F -(l-l3)eV) = 2iV_(e F ). 
This treatment is clearly not self-consistent and does not 
satisfy the condition of gauge invariance. We should also 
mention that we should treat the system in the grand 
canonical ensemble so that the Hartree potential on the 
force oscillation is not double-counted, as did in Ref. (23. 



IV. MAGNETIC PROPERTIES 

We have already seen that it is very important to 
treat the electron-electron interactions appropriately to 
calculate the cohesive force and the differential conduc- 
tance. The interactions can also produce significant ef- 
fects in magnetic properties. By using the free elec- 
tron model, the calculated magnetoconductance coeffi- 
cient a and magnetotension coefficient T in are identi- 
cally zeroi However, since the Hartree potential U has 
to be determined self-consistently by Eq. @ in presence 
of the magnetic field. As we will show that, the response 
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FIG. 3: (Color line) The magnetoconductance coefficient a 
versus the bias voltage. 
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FIG. 4: (Color line) The magnetotension coefficient T versus 
the bias voltage. 



of Hartree potential is sensitive to the external mag- 
netic field at thresholds of the single-particle subbands, 
and gives rise non-zero magnetoconductance coefficient 
and magnetotension coefficient at these sub-bands. We 
should also show that this response of Hartree potential 
to the magnetic field also affects the magnetic suscepti- 
bility. 

Using perturbation theory, the effect of a weak longitu- 
dinal magnetic field H (i.e. a field such that the cyclotron 
radius r c 3> Rq perpendicular to the cross section of the 
wire) can be included as a spin-dependent shift of the 
transverse eigenvalues^ 



gm s )— 

#0 



(9) 



where H = k F hc/2e, g is the gyro magnetic ratio factor, 
c is the speed of light, m v is the orbital angular momen- 
tum, and m s — ±^ is the spin of an electron. 

The generalized grand-canonical potential and the 
current at weak magnetic field are modified to be 



n [R 0l V,H] 



2e F L 
3A F 
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and 
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dE 0(E 



s). (11) 
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The Hartree potential U should be self-consistently de- 
termined by Eq. © under both voltage and the magnetic 
field. The cohesive force F and differential conductance 
G are also modified to be 



u.s 



e F 



1/2 



dRr 



and 
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5: (Color line) The derivative of the Hartree energy U 
respect to the magnetic field H versus the bias voltage. 



G{R 0l V,H) 



G 



E^ 

ot=± 
v.s 



dU 



(13) 

The calculated cohesive force and differential conduc- 
tance at small magnetic field from Eqs. i|12|) and i|13|) 
are not much different from Eqs. J7J and JSJ. These two 
equations serve the purpose to calculate the magnetoten- 
sion and magnetoconductance coefficients below. 

The magnetotension coefficient T is defined as T = 
(dF/dH)H=o and magnetoconductance coefficient a is 
defined as a — j^{dG/dH) H =o- Using Eqs. iTf^) and 
(IPf . one gets 
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6: (Color line) The magnetic susceptibility x versus volt- 
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part, one gets 



X (eV) 



E - 
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and 
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Xf 

a—± 
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SF 
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k 2 F Rl dL 



dU 
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H=0 



where g(e a ) is the density of states at energy e a of elec- 
trons injected from reservoir a. From these two equation, 
one can see that nonzero magnetotension and magneto- 
conductance coefficients are attributed to the response of 
the Hartree potential U to the magnetic field. 

The calculated results for a and T are shown in Figs. 
OH and 0] Whenever there is a subband opening, there 
is a peak in these two coefficients. This is because the 
magnetic field increases the fluctuation of the Hartree 
potential at those subbands, which can be seen in Fig. El 

One can consider the single-particle and interaction 
effects on the magnetic susceptibility which is define as 

X = — \ (IjJrt) h ■ Using the generalized grand canon- 
ical potential Eq. (|10fl . and keeping only the dominant 



d(U + e vs ) 
dH 



(16) 



The result is presented in Fig. For nanowires with 
small radii R , is of the same order as ^ff- = 
eF{m l/ +gm s )/Ho, both single-particle effect and interac- 
tion effects are important to the magnetic susceptibility. 
For nanowires with large Rq, Jjj- is negligible comparing 



to 



dH 



V. DISCUSSIONS AND CONCLUSION 

We should point out that based on our recent stabil- 
ity analysis^, the nonzero a and T and the spikes of G 
and x usually appear in the mechanical unstable zones. 
However, this mechanical instability should be a problem 
only in the measurement of T since such a measurement 
is a mechanical process. The sensitivity of the Hartree 
potential to the applied voltage and magnetic field can 
be observed by measuring the differential conductance G, 
magnetoconductance coefficients and magnetic suscepti- 
bility as a function of the applied voltage as long as the 
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measuring time is shorter than the lifetime caused by the 
mechanical instability. 

In conclusion, we have used the generalized mean field 
electron model at finite voltage biasi^ to analyze the 
single-particle and interaction effects on the cohesion, 
transport and magnetic properties of cylindrical metal 
nanowires. At finite voltage bias, it is crucial to treat 
the electron-electron interactions adequately so that the 
calculated physical quantities are gauge invariant. Our 
results show that the cohesive force is determined by the 
single-particle effect, while the nonzero magnetotension 
and magnetoconductance coefficients are attribute of the 



response of the Hartree potential to the magnetic field. 
Both single-particle and interaction effects are important 
to the differential conductance and the magnetic suscep- 
tibility. 
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